On the soiiton width in the incommensurate phase of spin-Peierls systems. 
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We study using bosonization techniques the effects of frustration due to competing interactions and 
of the interchain elastic couplings on the soiiton width and soiiton structure in spin-Peierls systems. 
We compare the predictions of this study with numerical results obtained by exact diagonalization 
of finite chains. We conclude that frustration produces in general a reduction of the soiiton width 
while the interchain elastic coupling increases it. We discuss these results in connection with recent 
measurements of the soiiton width in the incommensurate phase of CuGeOs. 
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One dimensional (ID) magnetic systems are unstable 
against a spin-lattice interaction. At a finite tempera- 
ture a phase transition takes place toward a dimerized or 
spin-Peierls (SP) state whose magnetic signature is the 
opening of a gap (SP gap) between the singlet grciund 
state and triplet excitations. The recent discoverytl of 
the first inorganic SP compound CuGeOa has renewed 
the interest in these quasi-one-dimensional spin-phonon 
coupled systems. A number of experiments can be per- 
formed on these inorganic materials in an easier way than 
on the organic SP systems. It is also now possible to re- 
place the Cu ions by other magnetic or nonraagnetic ions 
and to study the effects of these impurities.0 Adding in- 
terest to the study of this compound is the fact that, as it 
was early recognized, the temperature dependence of the 
spin susceptibility could only be accounted for if an im- 
portant next nearest-neighbor (NNN) is included in 
ID Heisenberg model used to describe this materiall 
The dispersion of magnetic excitations and a number of 
other experimentally determined properties [are also con- 
sistent with a model with NNN interaction.ULl 

Recent x-ray-scattering experiments have detected a 
transition from the uniform SP phase into an incommen- 
surate phasep in the presence of a magnetic field. The 
experimental x-ray spectra in the incommensurate state 
has been interpreted as due to a soiiton lattice structure. 
A soiiton deformation was early theoretically predicted 
for SP systemstl. It arises because the Zeeman energy fa- 
vors states with nonzero magnetic moment thus breaking 
the dimer pattern. That is, each down spin changed into 
an up spin breaks a singlet dimer and gives rise to two 
solitons, each soiiton carrying a spin-i. 

The soiiton width has been experimentally estimated 
from the relationship between the measwed intensities 
of the main x-ray peak and its harmonics.l3 The resulting 
value is ^ = 13.6c with c the lattice constant in the spin 
chain direction {c — 1 hereafter). The predicted width 
for a one-dimensional model with only nearest neighbor 
(NN) interaction is given by ^ = Jtt/ (2A), where J is the 
exchange coupling and A the SP gap.EI Using the experi- 
mental value of A = 2.1 meV and J — 120 K, estimated 
by fitting the spin susceptibility with a Heisenberg model 



with NN interactions, the predicted ^ is equal to 8 which 
is considerably smaller than the experimental one. 

It is then quite apparent that some other ingredients 
should be added to the model in order to explain the ori- 
gin of this disagreement between theory and experiment. 
One of the goals of the present work is to analize the ef- 
fect of the frustration induced by the competing NN and 
NNN antiferromagnetic interactions on the soiiton width. 
In this sense, we are going to extend the formalism which 
leads to the relationship between the soiiton width and 
the gap, originally developed in Ref. [ ^. 

Another possible contribution to the soiiton width 
comes from the interchain elastic coupling. It is well 
known that the three-dimensional (3D) character of the 
phonon field is essential to account for the finite tem- 
perature SP transition. Moreover, for the specific case 
of CuGeOa, it has been found that the principal lattice 
anomaly takes, place in a direction perpendicular to the 
spin chainBB and it is related with displacements of 
the oxygen ions toward and against the chains. In this 
situation the phonon modes perpendicular to the chains 
involving coherent motion of the atoms in different chains 
become relevant in the SP phase.EHl It is then reasonable 
to assume that the interchain elastic coupling could be 
relevant to explain some features of the system. The pos- 
sibility that interchain elastic coupling could modify the 
soiiton shape is supported by a recent measurement of 
the elastic constants in CuGeOa in the presence of high 
magnetic fields.E£l In this study, it was shown that the 
elastic constant along the b axis is the one that presents 
the largest variation at the dimerized-incommensurate 
transition. 

We conclude that frustration produces in general a re- 
duction of the soiiton width while the elastic coupling of 
a chain with its neighbors increases it. 

We start our study by analyzing the effects of frus- 
tration on the soiiton width. Therefore, we consider the 
following ID model Haniiltonian: 
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where Si are spin-1/2 operators and Ui is the displace- 
ment of magnetic ion i from its cquihbrium position. No- 
tice that in CuGeOs, the displacements Ui are related to 
the displacements of the O ions. Hei in Eq. (0) corre- 
sponds to the elastic energy interaction of the ions along 
the chain with dimensionless elastic constant K^\. We 
include the phonons in the adiabatic approximation and 
we allow arbitrary displacement patterns. 

We now go to the bosonic representation of Hamilto- 
nian (|^). This is achieved by first applying a Jordan- 
Wigner transformation on the spin operators and then 
linearizing the resulting fcrmionic relation dispersion 
around the Fermi level. Finally, the fcrmionic field theory 
is bosonized. Very precise predictions have been obtained 
within this approach for the ground state aniJpw energy 
excitations of one-dimensional spin systemsEJ. Nakano 
and FukuyamaH used this bosonic representation to study 
soliton formation in the unfrustrated chain. In the follow- 
ing, we extend their approach to the frustrated chain. In 
terms of the boson field (j) and its conjugate-momentum 
n the Hamiltonian (|l|) can be rewritten asO 



dxu 
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Hbos = 
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In going from (|^) to (||,^) we have replaced the lattice 
displacement Ui by a continuum field as: 



i-iyui u{x) 



(4) 



so that we can describe smooth deviations of the dis- 
placement field with respect to the dimerized pattern. 

The quantities Vg , t] and g are expressed in terms of the 
microscopic parameters of Hamiltonian (^. The expres- 
sions can be read in Ref. [ However, once identified 
the physical meaning of these quantities it is customapj; 
to replace those expressions by more precise relations. t2l 
In this sense, the isotropy of the long distance behav- 
ior of the correlation functions implies that ry = 2. The 
parameter g measures the strength of the umklapp pro- 
cesses in the fcrmionic theory. The last term of (^) is a 
marginally irrelevant operator for all g > 0, so that in 
absence of coupling to the lattice a gapless spin liquid 
phase is predicted. When g is negative the last term of 
Eq. (||) becomes a relevant operator leading-the system 
to a spontaneously dimerized gapped phasetJ Therefore, 



g should be proportional to ac — a where Ofc is the crit- 
ical value separating the gapless and the gapped phases. 
The quantity Vs is the spin wave velocity. For a = 0, 
Vs = Jtt/2 from Bethe's exact solution. For a < ac, 
has been recently numerically evaluated as the derivative 
of the-jnagnon dispersion relation at the bottom of the 
bandJl3 This derivative is not well defined for a > ac 
due to the presence of the gap above discussed. We will 
analyze later on our procedure to evaluate Vs in this case. 

The application of a magnetic field to the system fa- 
vors states with nonzero Sz and to the appearance of 
solitons. Each soliton could be regarded as a domain- 
wall separating two opposite dimerization patterns. To 
take into account these states we follow the approach of 
Ref. [ ||. We separate the variable </> into the sum of a 
classical variable 0o and its quantum fluctuation (j). A 
self-consistent harmonic approximation is then used. We 
retain (j) to quadratic order and we require the vanishing 
of the first order term. The displacement field u is ob- 
tained from its static equilibrium equation. The resulting 
equations are: 

- ^(9^<^) + u e-^-^'^/^ cos 00 - 25e-2<'^'> sin(2(/)o) = 



3in0oe"<'^'>/^ +4if||M = 
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where (0^) is the ground state expectation value which 
we assume the same as in the dimerized phase. These 
equations could be combined to give: 



The soliton solution of this equation is: 
sin((/)o) = tanh(x/^) 

The lattice soliton is then given by: 
u{x) = uotanh(a::/^) 



(6) 
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with Mo the equilibrium displacement in the uniform 
dimerized state. The main quantity we are interested 
in is the soliton width which is given by: 



1 

e 



V 

2v, 



-<ct>^>/2 



456 



(9) 



This quantity is related to the singlet-triplet gap in the 
uniform dimerized state in a simple way. In this state 
u{x) — uq and (po — 7r/2. By following similar steps 
as in Ref. [ ^ we obtain the dispersion relation for the 
excitations of the system given by: 



(10) 



where ko = l/S, is given by Eq. (g). Note that Vg is 
a function of a in this case. Finally, the simple scaling 
relation between the soliton width and the SP gap A, 
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(11) 



originally obtained for the unfrustrated chain, is also 
valid in the presence of frustration. For a > ac, fco 
contains a contribution from the frustration due to the 
presence of a gap even in the absence of dimerization. 

The relation (11) can be tested by numerical exact di- 
agonalization on finite chains. The adiabatic equations 
corresponding to Hamiltonian (|^) are: 



|S, -S.+ilO) =0 
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where Si ~ (wi+i — Ui) is the bond length variable. We 
have solved iteratively these equations for Si starting 
from a random configuration The ground state 

of the spin system |0} was recalculated at each iteration 
step by Lanczos diagonalization. We considered finite 
chains of even number of sites up to 20 sites with peri- 
odic boundary conditions. r£h& numerical details of the 
method are given elsewherell3. 

In the subspace = the system converges to a 
dimerized configuration {(^l"-*} as expected. For a = 0.0, 
0.2 and 0.4, and for several values of i^y , using the dimer- 
ized pattern obtained at — 0, we computed the spin 
gap as A = Eo{S^ = 1) - £^0(5"^ = 0). 

When the spin Hamiltonian is diagonalized in the sub- 
space 5*^ = 1 we found that two solitons appear in the 
dimer pattern. The deformation pattern {Si} could be 
described by a superposition of two solitons centered at 
sites ii and «2 of the form: 



5, = {-lySatanh 
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where Sq, the deformation amplitude, and ^, the soliton 
width, were obtained by numerical fitting. The main lim- 
itation of this calculation arises in the region where, for a 
given a, T^y is so large that the solitons have a substantial 
overlap in our small chain and the fitting function ( p^ ) 
is no longer appropriate. Unfortunately, it is precisely in 
this situation where the continuum approximation lead- 
ing to the relation ( |Tl| ) is expected to hold. 

In Fig. 1 we show the soliton width ^ as a function 
of 1/A for the same values of a as above. The linear 
behavior predicted by Eq. ( |ll]) is clearly seen. A lin- 
ear fitting of these curves in the region ^ > 2.5 gives 
the slopes 1.87, 1.70 and 1.63 , for a = 0.0, 0.2 aad. 
0.4 respectively. For a < Uc, a numerical calculation^ 
leads to Vg = ^(1 — 1.12q!) in the thermodynamic limit. 
From this law, Vg = 1.57, 1.22, for a = 0.0 and 0.2 re- 
spectively. We can observe that the slopes obtained by 
fitting the numerical data are larger than the calculated 
Vs- Besides, the effect of a is weaker in the numerical 
data than that predicted by Eq. (pj]). It is difficult 
to conclude at this point if this disagreement between 
the prediction obtained by the continuum bosonized the- 
ory and numerical results is due to the approximations 
involved in the former or to finite size effects. For 



a = 0.4 > ac ~ 0.2411, we have estimated Vg by fit- 
ting the energy e(fc)^ = A^ -I- v'^k'^ + ck^ around fc = 0, 
where s{k) = EqIs^ = l,fc) - £^0(5"^ = 0, fc = 0). For 
L = 20, we obtained Vs = 0.707, a value which is smaller 
than the slope of the curve (, vs. 1/A for a — 0.4 in Fig. 

®- . . . 

Then the qualitative picture of the effect of frustration 

is partially confirmed by the exact diagonalization study. 
Both analytically and numerically we have seen that frus- 
tration reduces the soliton width for a given value of 
the SP gap. However, the numerical calculations show 
a smaller dependence on a than that analytically pre- 
dicted. 
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FIG. 1. Soliton width vs. obtained by exact diagonal- 
ization in the 20 site chain for various values of a = J2/J. 

We now turn to the study of the effects caused by the 
3D character of the phonon field. As previously dis- 
cussed, there are both theoretical as well experimental 
arguments indicating that the interchain elastic coupling 
is important to describe the excitations of the system in 
the SP regime. Therefore, we add to H^i of Eq. (|l|) an 
interchain elastic coupling of the form: 



^ 2 ^ 



(14) 



where u\ denotes the displacement of the j ion in the 
chain j (j = {jx,jy))- K^^ and K'^ are the harmonic cou- 
plings in the a and b directions respectively. We assume 
that there is no interchain magnetic interactions. The 
system now consists of a set of spin chains immersed in 
a 3D lattice. We resort again to bosonization methods 
to analyze soliton formation in this system. The boson 
field (f> turns to be j-dependent. The adiabatic condition 
(^) is now given by: 

sin0j,e~<'^'>/2 + [2{Kl + K^) + AK\\]u^ - 
K l(u^-+i^^« - uJ--i'^») - Xl(uJ- - mJ-^«-1) = 0. 
This equation can be inverted to give: 
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^(^) = -I^EB(j'-j)sin0{(x)e- 
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where: 
B(j'-j) 



dk cos[k.(y-j)] 

^ 7-x (27r)2 1 + sin^ ^ + sin^ 



(16) 



where e"' = are the relative interchain elastic cou- 

plings. The equations governing the classical fields (jy^Q 
now become: 

5''/^o+2^B(O)sin(20j) 

+ l^B(j'-j)cos(/>jsin</;i' =0 (17) 

where we have only included for simplicity the effect of 
the NN exchange interaction in each chain. In this sit- 
uation, ^ is given by Eq. (^) with g = 0. Moreover, as 
the interchain elastic couplings do not affect the uniform 
dimerized state, the relation (^l|) between ^ and the gap 
is preserved. 

The second term of Eq. (|l^) is dominant due to the 
rapidly decaying effective interaction B(j' — j) with the 
distance |j' — j|. This behavior can be seen in Fig. Hfor 
the case e'^ — e^. By neglecting the last term in Eq. (|17|) 
we obtain a system of decoupled equations of the form 
(^. We consider a situation where a soliton is formed 
in a single chain in an otherwise SP state. Therefore we 
have 0j — Tr/2 for j ^ m and sm{(j>Q) — tanh(a;/fM for 
j = m. The lattice pattern is obtained from Eq. (15): 



= B(j - m)uo tanh(a;/C') + (1 - B(j - m))uo (18) 

This solution represents a soliton at the origin of the m 
chain. It is easy to see that the displacement amplitude 
at one side of the soliton is smaller than at the other side. 
This difference can be explained by the fact that at one 
side the dimerization has the same phase with respect to 
that of the neighboring chains while at the other these 
phases are opposite. At the same time, the neighbor 
chains j ^ m feel the presence of the soliton in m. As 
seen from Eq. ([l^ ) they are slightly distorted from the 
uniform dimerized pattern. 
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FIG. 2. The ratio between the effective soliton width and 
the "bare" soliton width as a function of e" = e*" (circles). 
Coefficients B(0) (filled triangles) and B{1,0) (open triangles) 
given by Eq. (|lq) are also shown. 

From E q. (p^ ) we obtain the effective soliton width 
C = ^/ -^5(0) which is larger than the "bare" width 
^ of the one chain problem previously discussed since 
B{0) < 1. We show in Fig. || the behavior of ^'/^ as 



a function of e" for the case of e" 



It is difficult 



to extract from experimental data an estimation of the 
values of e°"^ since these are effective elastic constants as 
introduced in the interaction term (p^). 

After this work has been-eompleted we became aware 
of a paper of Zang et alSJi where the effects of frus- 
tration and interchain magnetic coupling on the soliton 
width were studied. Their conclusion about the frustra- 
tion effects is similar to the one presented here. Besides, 
they conclude that the large soliton width experimentally 
found in CuGeOa is mainly due to the antiferromagnetic 
interchain coupling. However, as the study of Ref. [ |l^] 
is based on a mean-field approximation for the interchain 
coupling and since the soliton formation is a strongly lo- 
cal phenomena, we think that the enhancement of the 
soliton width found by these authors might be overesti- 
mated. In this contribution from the interchain 
elastic coupling could be needed to explain the experi- 
mental result. The question of to what extent these two 
mechanisms contribute to the soliton width could only 
be answered by studying a model which includes both 
magnetic and elastic interchain couplings. 
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